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ABSTRACT 

We observed two fields near M32 with the Advanced Camera for Surveys/High Resolution Channel 
(ACS/HRC) on board the Hubble Space Telescope (HST). The main field, Fl, is 1.'8 from the center of M32; 
the second field, F2, constrains the M31 background, and is 5. '4 distant. Each field was observed for 16-orbits 
in each of the F435W (narrow B) and F555W (narrow V) filters. The duration of the observations allowed 
RR Lyrae stars and other short-period variables to be detected. A population of RR Lyrae stars determined to 
belong to M32 would prove the existence of an ancient population in that galaxy, a subject of some debate. 

We detected 17 RR Lyrae variables in Fl and 14 in F2. A Id upper limit of 6 RR Lyrae variables belonging to 
M32 is inferred from these two fields alone. Use of our two ACSAVFC parallel fields provides better constraints 
on the M31 background, however, and implies that 7^ (68 % confidence interval) RR Lyrae variables in 
Fl belong to M32. We have therefore found evidence for an ancient population in M32. It seems to be 
nearly indistinguishable from the ancient population of M31. The RR Lyrae stars in the Fl and F2 fields have 
indistinguishable mean V-band magnitudes, mean periods, distributions in the Bailey diagram and ratios of 
RRc to RRtotai types. However, the color distributions in the two fields are different, with a population of red 
RRab variables in Fl not seen in F2. We suggest that these might be identified with the detected M32 RR Lyrae 
population, but the small number of stars rules out a definitive claim. 

Subject headings: Local Group — galaxies: individual: M32, M31 — galaxies: elliptical and lenticular, cD — 
stars: Population II — stars: variables: other 



1. INTRODUCTION 

Messier 32 (M32) is the only elliptical galaxy close enough 
to possibly allow direct observation of its stars down to the 
main-sequence turn-off (MSTO). It is a vital laboratory for 
deciphering the stellar populations of all other elliptical galax- 
ies, which can only be studied by the spectra of their inte- 
grated light, given their greater distances. Major questions 
about M32's star formation history remain unanswered. M32 
appears to have had one or more relatively rece nt episodes 
of star formation (with i n the last 3 Gyr: e.g.. lO'Connelll 



oi star tormation (witmn tne last i Uyr: e.g.. IU C onnell 
19801: lRosdll98lJ[l994t lGonzaledll993t iTrager et al.ll200oT 
Coelh o et al.l 2009), whi ch also appears t o be true for many 
elliptical galaxies (e.g., lGonzaledll993t ITrager etafl l2000t 

Electronic address: fiorentino@astro.rug.nl 

1 Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations are associated with GO 
proposal 10572. 

2 The National Optical Astronomy Observatory is operated by AURA, 
Inc., under cooperative agreement with the National Science Foundation. 



iThomas et alj 12005). These conclusions rest on painstak- 
ing and controversial spectral analysis of their integrated 
light. In contrast, the most direct information about a stel- 
lar population comes from applying stellar evolution theory to 
color-magnitude diagrams (CMDs). Little however is known 
about M32's ancie nt population (see, e.g., iBrown et alj|2000t 
iCoelho et al.ll20 09). 

With our Advanced Camera for Surveys/High Resolution 
Channel (ACS/HRC) data (Cycle 14, Program GO- 10572, 
PI: T. Lauer) we have obtained the deepest CMD of M32 to 
date. A compreh ensive analysis of this C MD is discussed in a 
companion paper Monac hesi et al.l (120091 hereafter M09) and 
we refer to it for further details. However here we want to 
stress that, due to the severe crowding in our fields, even with 
the high spatial resolution of HRC it is not possible to reach 
the MSTO with sufficient precision to claim the presence of a 
very old population. 

RR Lyrae variables are low-mass stars burning He in their 
cores. They are excellent tracers of ancient stellar popula- 
tions, completely independent of the MSTO, and knowledge 
of their properties provides important information on their 
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parent stellar populations. Because they are located on the 
horizontal branch (HB) in a CMD, they are at least 3 magni- 
tudes brighter than MSTO dwarfs and therefore detectable to 
relatively large distances. RR Lyrae are also very easy to char- 
acterize, with ab-type RR Lyrae (RRab) pulsating in the fun- 
damental mode, rising rapidly to maximum light and slowly 
declining to minimum light, and c-type RR Lyrae (RRc) pul- 
sating in the first harmonic mode, with their luminosities vary- 
ing roughly sinusoidally. Most importantly for our purpose, 
the mere presence of RR Lyrae stars among a population of 
stars suggests an ancient origin, as ages older than ~ 10 Gyr 
are required to produce RR Lyrae variables. Thus the detec- 
tion of RR Lyrae stars in M32 is presently the only way to 
confirm the existence of an ancient stellar population in this 

g alaxy. 

lAlonso-Garcfa et al] J2004) were the first to attempt to di- 
rectly detect RR Lyrae stars in fields near M32. They imaged 
a field ~ 3 '5 with WFPC2 to the east of M32 and compared 
it with a control field well away from M32 that should sam- 
ple the M31 field stars. They identified 12 ± 8 variable stars 
claimed to be RR Lyrae stars belonging to M32 and therefore 
suggested that M32 possesses a population that is older than 
~ 10 Gyr. They were however unable to classify these RR 
Lyrae variables and could not derive periods and amplitudes 
for them. 

Very recently, Saraiedi ni et al.l ([2009, hereafter S09) used 
ACS/WFC parallel imaging from our present data set to find 
RR Lyrae variables in two fields close to M32. They found 
68 1 RR Lyrae variables, with excellent photometric and tem- 
poral completeness (Sec. These RR Lyrae stars were 
roughly equally distributed between the two fields, with the 
same mean average magnitudes, metallicities, and Oosterhoff 
types in each field. It was therefore impossible for them to 
separate the variables into M31 and M32 populations. It is 
still therefore an open question as to the precise nature or 
even presence ofRR Lyrae variables in M32. 

In this paper, we present newly-detected RR Lyrae vari- 
ables observed with ACS/HRC and also a detailed analysis 
of the fields near M32 where RR Lyrae stars have been found 
with HST The paper is organized as follows. In Section 2 
we describe our observations and the data reduction we per- 
formed. We move on to describe the technique used to iden- 
tify and characterize the RR Lyrae variable stars in Section 3, 
where we present their periods and light curves. In Section 4 
we show that we have clearly detected RR Lyrae variables in 
M32, as long as we include the results from our ACS/WFC 
parallel fields. In Section 5 we discuss the properties of the 
RR Lyrae stars, such as the location of their instability strip, 
reddenings, mean periods and Oosterhoff types, as well as 
pulsational relations such as period-metallicity-amplitude. In 
this section we also derive estimates of the distance moduli to 
and metallicities of our fields. We summarize our findings and 
present our final conclusions in Section 6. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Field Selection, Observational Strategy and Data 
Reduction 

We obtained deep B and V-band imaging of two fields near 
M32 using the ACS/HRC instrument on board HST during 
Cycle 14 (Program GO-10572, PI: Lauer). The primary goal 
of this program was to resolve the M32 MSTO. The ACS 
F435W (B) and F555W (V) filters were selected to optimize 
detection of MSTO stars over the redder and more luminous 



stars of the giant branch. M32 is very compact and is pro- 
jected against the M3 1 disk. Thus the major challenge was 
to select a field that represented the best compromise between 
the extreme crowding in M32, which would drive the field to 
be placed as far away from the center of the galaxy as possi- 
ble, versus maximizing the contrast of M32 against the M3 1 
background populations, which would push the field back to- 
wards the central, bright portions of M32. Following these 
constraints, the M32 HRC field (designated Fl) was centered 
on a location 1 10" south (the anti-M31 direction) of the M32 
nucleus, roughly on the major axis of the galaxy. The V- 
band surface brightness of M32 ne ar the center of the field 
is fly w 21.9 (Kor mendv et al.ll2009l) . M32 quickly becomes 
too crowded to resolve faint stars at radii closer to the center, 
while the galaxy rapidly falls below the M3 1 background at 
larger radii. 

Even at the location of Fl, M31 contributes ~ 1/3 of the 
total light, thus it was critical to obtain a background field, 
F2, at the same isophotal level in M31 (pv ~ 22.7) to allow 
for the strong M3 1 contamination to be subtracted from the 
analysis of the M32 stellar population. F2 was located 327" 
from the M32 nucleus at position angle 6 5° . At this angula r 
distance M32 has an ellipticity e ps 0.25 (IChoi et all 2002). 
and F2 is nearly aligned with the M32 minor-axis. Thus 
the implied semi-major axis of the M32 isophote that passes 
though F2 is 435", significantly larger than the nominal angu- 
lar separation. The estimated M32 surface brightness at F2 is 
Pv w 27.5, based o n a modest extrap olation of the B-band sur- 
face photometry o flChoi et alJ (12002b and an assumed color of 
B-V w 0.9. The contribution of M32 to F2 thus falls by a fac- 
tor of ~ 1 80 relative to its surface brightness at Fl . While one 
might have been tempted to move F2 even further away from 
Fl, it clearly serves as an adequate background at the location 
selected, while uncertainties in the M31 background would 
increase at larger angular offsets. The locations of both the 
Fl and F2 fields are shown in Figure [T] 

Detection of the MSTO required deep exposures at Fl. Ac- 
curate treatment of the background required equally deep ex- 
posures to be obtained in F2. A summary of the observations 
is shown in Table [2 briefly, each field was observed for 16 
orbits in each of the F435W and F555W filters for a total 
program of 64 orbits. While the detection of RR Lyrae vari- 
ables was not the primary goal of the program, execution of 
each filter/field combination in a contiguous time span of 2- 
3 days was clearly well-suited to detect RR Lyrae variables, 
which have periods ranging from 0.2-1 d. 

At B and V, the HRC undersamples the PSF, despite its 
exceptionally fine pixel scale. All of the images were ob- 
tained in a 0.5 x 0.5 sub-pixel square dither pattern to ob- 
tain Nyquist sampling in the complete data set. In detail, the 
sub-pixel dither pattern was executed across each pair of or- 
bits, with each orbit split into two sub-exposures. The tele- 
scope was then offset by 0."125 steps between the orbit pairs 
in a "square-spiral" dither pattern to minimize the effects of 
"hot pixels," bad columns, and any other fixed-defects in the 
CCD, on the photometry at any location. The data for each fil- 
ter/field combination thus comprises 8 slightly different point- 
ings, with Nyquist-sampling obtained at each location. In 
practice, the dithers were extremely accurate, and Nyqu ist im- 
ages c ould readily be constructed using the algorithm o fiLauerl 
(U9991) . 

Our optimal average photometry has been obtained 
by using these very deep, super-resolved images. We 
have performed photometry by using both the DAOPHOT 
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TABLE 1 
Log of observations 



Field 


O/2000.0 


5/2000.0 


Filter 


Exposure time 


Date 


Fl 


00 42 47.63 


+40 50 27.4 


F435W 


16x1279+ 16x1320 


20-22 Sep 2005 


Fl 


00 42 47.63 


+40 50 27.4 


F555W 


16x1279+ 16x1320 


22-24 Sep 2005 


F2 


00 43 7.89 


+40 54 14.5 


F435W 


16x1279 + 16x1320 


6-8 Feb 2006 


F2 


00 43 7.89 


+40 54 14.5 


F555W 


16x1279+ 16x1320 


9-12 Feb 2006 



K r • 




FIG. 1 . — Location of our pointings near M32 observed with ACS/HRC on board HST (small squares). Fields Fl and F2 are at distances of 1 f 8 and 5'4 from 
the center of M32, respectively, and each of them covers a region of 0.25 arcmin 2 on the sky. Parallel fields, taken with ACS/WFC, F3 and F4 are at distances of 
5'3 and 9'2 from the center of M32, respectively, covering a region of 9 arcmin 2 on the sky. 



II/ALLSTAR packages dStetsonlll987l [T994) and by first de- 
convolving those combined images with a reliable PSF and 
then performing aperture photometry on the deconvolved im- 
ages. Both methods returned comparable results and allow us 
to present the deepest CMD of M32 obtained so far. This re- 
sult is analyzed in M09 and will not be discussed further here. 
However, in what follows we will use these CMDs to show 
the location of RR Lyrae stars. 

In addition to the HRC images, parallel observations were 
obtained with the ACS/WFC channel using the F606W filter 
(broad V). These fields, designated F3 and F4, are also shown 
in FigureQ] Notably, the telescope rolled by roughly 180° be- 
tween the execution of the Fl and F2 observations, thus F3, 
the parallel field associated with Fl, and F4, the mate to F2, 
bracket the Fl and F2 fields in angle. By happenstance, the 



F3 and F4 fields also nearly fall on the same M31 isophote 
that encompasses the Fl and F2 fields, thus the M31 back- 
ground should be roughly similar in all four fields. It is also 
notable that F3 is positioned slightly closer to the M32 nu- 
cleus than F2 (317" versus 327"), but because it also falls 
along the M32 major rather than minor axis, its associated 
M32 surface brightness is /jy ~ 25, or a factor of ~ lOx more 
than the M32 contribution to F2. Furthermore we note that the 
parallel observations were exposed only at the same time as 
the F555W exposures in Fl and F2 and therefore cover only 
half of the total time window of the primary exposures (we 
return to this point in Sec. |4]below). 

The parallel images have already been analyzed by S09. 
They find 68 1 RR Lyrae variables stars, of which 324 are lo- 
cated in the field closest to M32. Because only one filter was 
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available for the parallel observations, S09 did not have all the 
information needed to properly disentangle the populations 
that belong to M31 and/or M32. In fact, their detected RR 
Lyrae stars show the same mean average magnitude, metal- 
licity and Oosterhoff type, as we discuss in Section [5] With 
our primary observations we can attempt to disentangle the 
two populations by using all the quantities characterizing the 
class of RR Lyrae variables, such as mean weighted magni- 
tudes and colors in the Johnson-Cousins system, periods, and 
amplitudes. 

2.2. Photometry of the RR Lyrae Variables 

The study of the presence of RR Lyrae stars is based on a 
detailed analysis of the time series of our fields. We analyzed 
each single epoch image (32 per field and per filter) and not 
the combination of all the images described above. Because 
of the intrinsic brightness of the RR Lyrae (V ~ 25 mag), we 
decided to perform PSF-fitting photometry over all the fully 
calibrated data products (FLT) i mages using th e DOLPHOT 
package, a version of HSTphot (Dolphin 2000) modified for 
ACS images. Our choice has been justified by the short time 
consumed to obtain high quality photometry at the RR Lyrae 
magnitude level for our data set. Following the DOLPHOT 
User's Guide, we have performed the pre-processing steps 
MASK and CALCSKY routines before running DOLPHOT. 
This package performs photometry simultaneously over all 64 
images of each field, returning a catalog of more than 20000 
stars per field already corrected for charge transport efficiency 
(CTE) and aperture correction, following the suggestions by 
Siria nni et all (120051) . We use this photometry to perform the 
analysis of variable stars. 

2.2.1. Photometric Completeness at the Horizontal Branch 

The 50% completeness levels of fields Fl and F2 are at least 
2 mag deeper than the horizontal branch, as shown by the 
artificial star tests (ASTs) in M09. The ASTs show that the 
completeness at the red clump (i.e., the red horizontal branch) 
is 100%. The ASTs in M09 do not properly populate the blue 
horizontal branch, as there are very few stars in this region of 
the CMD, so we assume that the completeness at the position 
of the blue horizontal branch is the same as at the red clump. 

3. SEARCHING FOR VARIABLE STARS AND THEIR PERIODS 

To identify variable sources in both fields, we used a 
code written by one of us (AS) in the Interactive Data Lan- 
guage (IDL) whose pri nciples, based on t he algorithm of 
iLafler & Kinmanl (119651) . are discussed in Sah a & Hoessell 
( 1990). This code was applied to the results from DOLPHOT 
PSF-fitting photometry described in the previous section. Its 
output gives us not only a list of candidate variable stars but 
also a good initial estimate of their periods. The method as- 
sumes that realistic error estimates for each object at each 
epoch are available from the photometry, which are first used 
to estimate a chi-square based probability that any given ob- 
ject is a variable. A list of candidates is then chosen, and each 
candidate is tested for periodicity and plausible light curves. 
The graphical interface of this program clearly shows possi- 
ble aliases and allows the user to examine the light curves 
implied for each such alias. The final decision making is done 
by the user. A refinement of the period for all the candidate 
variables has been performed by using two other independent 
codes. We used the Period Dispersion Minimization (PDM) 
algorithm i n the IRAF enviro nment to confirm the found pe- 
riodicity ( S telling wern 1 1 9781) . Further refinement was then 



obtained by using GRATIS (GRaphical Analyzer of Time Se- 
rie s, developed by P. Mont egriffo at the Bologna Observatory; 
see lClementini et aT] |2000 and references therein for details), 
which permits us to fit Fourier series to the magnitudes in each 
passband as a function of their phase. 

The magnitudes returned by DOLPHOT have already been 
calibrated onto the HST VEGAMAG photometric system, but 
for the following analysis we need to transform them onto the 
Johnson-Cousins (JC) system. We need therefore to take into 
account the color variations in the periodic cycles of the vari- 
able stars. We thus associate each phased epoch in the F435W 
filter with the corresponding best-fitting F555W model pro- 
vided by GRATIS at th a t epoc h, and vice versa. Finally, we 
apply the lSirianni et al.l d2005) transformations from F4-35W 
and F555W to B and V for ACS/HRC, and we re-analyze the 
new JC time series with GRATIS to improve the light curve 
models as well as the previously-constrained periods. The or- 
der of the Fourier series used to obtain the best fit and the 
epoch corresponding to maximum of the light curve at phase 
0=1 are given in columns 5 and 6 in Tables [2] and [3] This 
procedure allows us to derive well-sampled and consistent 
light curves in both filters. Proper periods, mean magnitudes 
weighted on the proper light curve in both HST VEGAMAG 
and JC photometric systems, colors, and amplitudes are given 
in Tables |2]and|3]for a total number of 31 bona fide RR Lyrae 
stars: 17 in Fl and 14 in F2. The time series photometry is 
given in Table [4] 

Finding charts for the newly-detected RR Lyrae variables 
are shown in Figures [2] and [3] on the combined F555W im- 
ages of Fl and F2, respectively. Their locations are shown 
in the CMDs from M09 calibrated onto the HST VEGAMAG 
photometric system using the average magnitudes as reported 
in Tables [2] and [3] (Figure [4j. We cross-correlated the RR 
Lyrae coordinates and magnitudes with the photometric cata- 
log from M09 to confirm the presence of our new RR Lyrae 
stars in the average photometry. For all the RR Lyrae variables 
we found a star with same coordinates and similar magnitude 
(see the zoomed-in CMDs in Fig.|5]l. 

An atlas of the light curves for all the newly-detected RR 
Lyrae variables is shown in Figures [6] and [7] In the atlas 
both data points in Johnson-Cousins system as well as the 
models used to perform a proper calibration onto this pho- 
tometric system are shown. The error bars take into account 
both the scatter between the data and the model used to fit the 
Fourier series (see columns 13 and 14 in Tables [2] and [3) and 
the photometric errors as returned by DOLPHOT program. 
By averaging the Johnson-Cousins magnitudes we have ob- 
tained (V) =25.34±0.15magforFl and (V) = 25.30±0.12 
mag for F2. Then, we have classified RR Lyrae variables into 
fundamental-mode (FU) or first-overtone (FO) pulsators by 
an inspection of this atlas. FO pulsators have mean periods of 
~ 0.3 d and sinusoidal light curves, whereas the FU pulsators 
have longer periods ( (P a b) = 0.59 ±0.11 d) and more com- 
plicated light curves (with up to 4 harmonics). Because our 
sample of RR Lyrae variables is small, we find the same ratios 
of FO to FU pulsators in the two fields to within the Poisson 
errors: N c /N toM = 0.23!°;g and N c /N toM = 0.42^;| for Fl 
and F2, respectively 3 . We discuss the RR Lyrae properties in 

3 The central estimates and la confidence intervals on N c /N\oal have been 
found using simulations of 10 5 Poissonian deviates of N totil \ and N c in each 
case. We use the median of the resulting distribution of the ratio values as the 
central estimate and the region of the diagram that contains 68% of the area of 
the probability distribution function to compute the la confidence intervals. 
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TABLE 2 
Fl RR Lyrae properties 



Star id 


RA 


DEC 


Period 


H a 


Epoch^ 


F435W 


F555W 


(V) 

\ f I 


(B)-(V) 


A B 


Ay 


resB c 


resV c 


tvDe 




J2000 


J2000 


(days) 




(JD) 


mag 


mas 


mag 


mag 


mag 


mag 


mag 


mag 




i 


00:42:48.435 


+40:50:32.11 


0.255 


1 


2453634.900 


25.56 


25.39 


25.44 


0.19 


0.64 


0.59 


0.11 


0.10 


RRc 


2 


00:42:46.509 


+40:50:30.19 


0.285 


1 


2453633.060 


25.68 


25.45 


25.50 


0.24 


0.77 


0.56 


0.12 


0.11 


RRc 


3 


00:42:46.563 


+40:50:23.03 


0.311 


1 


2453632.810 


25.42 


25.22 


25.30 


0.19 


0.61 


0.58 


0.09 


0.11 


RRc 


4 


00:42:48.237 


+40:50:12.85 


0.317 


1 


2453632.900 


25.59 


25.43 


25.47 


0.20 


0.61 


0.42 


0.12 


0.10 


RRc 


5 


00:42:48.384 


+40:50:30.61 


0.475 


2 


2453636.290 


25.52 


25.29 


25.34 


0.25 


0.99 


0.85 


0.10 


0.11 


RRab 


6 


00:42:47.349 


+40:50:41.91 


0.486 


3 


2453634.450 


25.93 


25.56 


25.60 


0.41 


1.09 


0.82 


0.15 


0.15 


RRab 


7 


00:42:48.074 


+40:50:30.17 


0.519 


4 


2453634.415 


25.34 


25.04 


25.09 


0.31 


1.01 


0.75 


0.08 


0.06 


RRab 


8 


00:42:47.552 


+40:50:30.41 


0.521 


3 


2453634.560 


25.68 


25.41 


25.46 


0.28 


0.93 


0.76 


0.09 


0.10 


RRab 


9 


00:42:47.087 


+40:50:43.08 


0.523 


3 


2453635.650 


25.61 


25.25 


25.29 


0.40 


1.09 


1.06 


0.12 


0.13 


RRab 


10 


00:42:47.034 


+40:50:28.29 


0.546 


3 


2453632.216 


25.93 


25.51 


25.55 


0.45 


1.41 


1.03 


0.18 


0.16 


RRab 


11 


00:42:47.462 


+40:50:42.81 


0.564 


2 


2453635.490 


25.51 


25.19 


25.21 


0.36 


0.74 


0.68 


0.15 


0.09 


RRab 


12 


00:42:48.737 


+40:50:32.31 


0.621 


2 


2453632.600 


25.43 


25.19 


25.18 


0.30 


1.25 


0.89 


0.17 


0.15 


RRab 


13 


00:42:47.014 


+40:50:36.27 


0.625 


2 


2453637.280 


25.57 


25.28 


25.29 


0.35 


0.88 


0.68 


0.09 


0.08 


RRab 


14 


00:42:46.411 


+40:50:29.96 


0.626 


3 


2453632.500 


25.53 


25.15 


25.19 


0.43 


1.09 


0.81 


0.09 


0.09 


RRab 


15 


00:42:47.554 


+40:50:16.48 


0.645 


3 


2453634.680 


25.83 


25.47 


25.50 


0.40 


0.74 


0.39 


0.10 


0.07 


RRab 


16 


00:42:46.075 


+40:50:25.38 


0.728 


2 


2453637.320 


25.74 


25.22 


25.25 


0.58 


0.39 


0.38 


0.12 


0.08 


RRab 


17 


00:42:46.638 


+40:50:25.24 


0.851 


2 


2453637.693 


25.60 


25.15 


25.17 


0.49 


0.65 


0.45 


0.10 


0.08 


RRab 



d Order of the Fourier series used to obtain the best fit 

k Julian Date where each curve shows its maximum of light at phase (p = 1 

c RMS deviation of the data points from the fitting model, in B- (res5) and V- (resV 7 ) bands respectively 



TABLE 3 
F2 RR Lyrae properties 



Star id 


RA 


DEC 




Period 


H a 


Epoch b 


F435W 


F555W 


(V) 


(B)-(V) 


A B 


A v 


resB c 


resV c 


type 




J2000 


J2000 




(days) 




(JD) 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 




1 


00:43:07.766 


+40:54:15. 


.31 


0.267 


2 


2453772.400 


25.46 


25.32 


25.38 


0.20 


0.66 


0.59 


0.09 


0.09 


RRc 


2 


00:43:07.704 


+40:54:23. 


18 


0.287 


2 


2453774.700 


25.47 


25.36 


25.41 


0.18 


0.70 


0.55 


0.10 


0.07 


RRc 


3 


00:43:08.518 


+40:54:08. 


86 


0.320 


2 


2453774.125 


25.51 


25.26 


25.31 


0.33 


0.46 


0.43 


0.08 


0.09 


RRc 


4 


00:43:07.734 


+40:54:27. 


17 


0.326 


2 


2453774.310 


25.51 


25.36 


25.41 


0.24 


0.57 


0.47 


0.10 


0.09 


RRc 


5 


00:43:07.874 


+40:54:31. 


99 


0.350 


1 


2453777.220 


25.45 


25.25 


25.30 


0.28 


0.60 


0.42 


0.11 


0.09 


RRc 


6 


00:43:08.300 


+40:54:21. 


88 


0.383 


2 


2453771.830 


25.25 


25.08 


25.14 


0.23 


0.63 


0.49 


0.07 


0.07 


RRc 


7 


00:43:07.671 


+40:54:01. 


99 


0.482 


3 


2453778.050 


25.53 


25.33 


25.39 


0.28 


1.37 


1.08 


0.15 


0.12 


RRab 


8 


00:43:08.339 


+40:54:23. 


36 


0.502 


3 


2453774.324 


25.64 


25.32 


25.37 


0.39 


1.19 


1.05 


0.12 


0.14 


RRab 


9 


00:43:07.594 


+40:54:21. 


32 


0.528 


3 


2453773.198 


25.53 


25.27 


25.30 


0.35 


1.27 


1.02 


0.11 


0.11 


RRab 


10 


00:43:07.424 


+40:54:26. 


69 


0.528 


3 


2453774.846 


25.40 


25.17 


25.23 


0.29 


1.26 


0.89 


0.13 


0.12 


RRab 


11 


00:43:07.727 


+40:54:13. 


25 


0.571 


2 


2453775.520 


25.52 


25.28 


25.35 


0.28 


1.06 


0.99 


0.15 


0.14 


RRab 


12 


00:43:08.107 


+40:54:21. 


00 


0.588 


3 


2453775.045 


25.64 


25.36 


25.41 


0.35 


0.94 


0.81 


0.10 


0.10 


RRab 


13 


00:43:07.977 


+40:54:21. 


33 


0.697 


2 


2453778.053 


25.40 


25.12 


25.19 


0.33 


0.70 


0.58 


0.09 


0.06 


RRab 


14 


00:43:07.798 


+40:54:14. 


39 


0.790 


2 


2453778.450 


25.28 


24.95 


25.00 


0.40 


0.63 


0.56 


0.07 


0.08 


RRab 



a Ol der of the Fourier series used to obtain the best fit 

k Julian Date where each curve shows its maximum of light at phase 0=1 

c RMS deviation of the data points from the fitting model, in B- (resS) and V- (resV 7 ) bs 

detail in Section[5] 

Here we want to conclude by addressing an important ques- 
tion about the temporal completeness of these observations. 
That is, could we have detected all of the RR Lyrae stars 
in these fields at any reasonable period? To compute the 
probability of detecting variability wit h periods of 0.2- 2 d, 
we followed the metho d suggested by ISaha et all d 19861) and 
ISaha & Hoessell d 19901 see section IV and Figure 7 therein), 
using software kindly supplied by E. Bernard. We simulated 
one million stars randomly phased and distributed with pe- 
riods of 0.2-2 days in bins of 0.001 d and then folded the 
Heliocentric Julian Dates of both filter datasets according to 
the random period a nd initial phase of each artificial star (see 
iBernard eT al. 2009, for details). A variable is considered re- 
covered if it has a) at least two observations around the max- 
imum of the light curve, b) at least two phase points in the 
descending part of the light curve, and c) a minimum of three 
observations during the minimum light. The results are shown 



s respectively 

in Figure [8] where the probability of detecting variability is 
plotted as a function of the period. The computed probabil- 
ity for Fl and F2 is nearly unity over the entire range 0.2-0.9 
days. We therefore assume a final (photometric plus tempo- 
ral) completeness parameter of 100% for both fields. 

4. HAVE WE DETECTED M32 RR LYRAE VARIABLE STARS? 

We have clearly detected RR Lyrae variables in Fl, a field 
dominated by M32. Are any of these stars truly associated 
with M32? Or does the strong background signal from M31 
RR Lyrae variables (judging from F2, which is nearly free of 
M32 stars) dominate our detection? 

We begin addressing these questions by examining the im- 
plications of our detections of RR Lyrae variables in Fl and 
F2 on the detection of M32 RR Lyrae variables. We then 
extend our analysis to include the M31 background repre- 
sented by fields F2-F4 and ask this question again. Fi- 
nally, we examine our re sults in the context of the study of 
lAlonso-Garcia et all (120041) . who have previously claimed de- 
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FIG. 2. — Finding chart for RR Lyrae variables found in field Fl. Numbers correspond to variables listed in Table|2] The FoV size is 0.25 arcmin 2 , as listed 
in column 9 of Table[5] We note that the R R Lyr ae stars in Fl are slightly clustered along the edge closest to the center of M32 where the total stellar density is 
increasing; we return to this point in Section l5~ll The arrow points towards the center of M32. 



tection of M32 RR Lyrae variable stars and therefore the pres- 
ence of an ancient stellar population in that galaxy. 

4.1. M32 RR Lyrae Population Inferred from Fl and F2 

We ask the question whether any of the RR Lyrae stars in Fl 
could belong to M32 given the M3 1 background represented 
by F2, or at least what the upper limit on the number of RR 
Lyrae stars belonging to M32 is. We have observed 14 RR 
Lyrae stars in F2 and 17 in Fl. Our assumption is that the 
stellar population in F2 represents a constant background in 
Fl. Let us call the "true" number (the expectation value) of 
RR Lyrae stars in each field belonging to M31 fig and the 
"true" number of RR Lyrae stars in Fl belonging to M32 fi F . 
What we observe is N B = 14 in F2 and Nb+N f = 17 in Fl . 
Then, by Bayes' theorem (see, e.g. ISivia & Skilling 2006), 
the probability of finding some /i# given N B is 4 



P(Vb\N b )cxP(Nb\»b)P(»b), 



(1) 



4 Note that here we are surpressing the role of the background information 
/, so that, for example, P(Nb) is shorthand for P(Nb \I). 



where the constant of proportionality, 1 /P{Nb), can be treated 
as a normalization constant such that J P(nB\Ns)d /ig = 1. By 
the product rule, the joint probability of finding fi F and /!# 
given Nb+Nf is 

P(^f^b\Nb+N f )^P(Nb+N f \p f ,pb)P(m\Nb)P(Pf), (2) 

where P(/j,b) and P(/j, f ) are priors on /j,b and p F which we 
discuss below, and P{Nb\ij.b) and P{NB + N F \/ji F ,fiB) are both 
represented by Poisson distributions, as we are counting stars: 



P(N\n) : 



H N e~v 
AH 



(3) 



In order to determine the probability distribution of /j, F given 
N f + Nb, we need to marginalize equation §2^ over /ig: 

P(p F \N B +N F )K J dfi B P(NB+N F \n F .jiB)P(^B\N B )P(fi F ). 

(4) 

All that is left is to specify the priors P(hb) and P(/j, f ) (and 
to normalize the probability distributions). Because this is a 
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FIG. 3. — Findingcharts for RR Lyrae variables found in field F2. Numbers correspond to variables listed in Table|3] The FoV size is 0.25 arcmin 2 , as listed in 
column 9 of Table|5] The adopted intensity scale is the same as in Fig.[2]for a fair comparison. 



location problem (see Chap. 5 of iSivia & Skilling 20061), we 
choose uniform priors, with ranges specified by the reason- 
able range of RR Lyrae specific frequencies Srr given the 
known old populations in M32 and M3 1 from M09 [see equa- 
tion[8]below]: 

W \ otherwise ' K } 

where a and b define the range of reasonable values for the 
expectation value /x. For fip, we choose a = and b = 25, 
and for ps, we choose a = 1 to enforce the presence of some 
background stars in Fl and b = 25. These limits on the pri- 
ors imply < S RR < 23 for M32 in Fl and 1.8 < S RR < 44 
for M3 1 (assuming the M3 1 stellar population is identical in 
Fl and F2) using the old, metal-poor star fractions described 
below. Note that we require integer numbers of stars, so the 
integration over \i% in equation © becomes a sum over /xg, 
restricted by equation (0 and our choice of a and b to the 
range 1 < p B < 25. 



With this machinery in place, we find the most probable 
[if = 3 RR Lyrae stars (as expected) belonging to M32, with 
68% confidence interval of 0-6 RR Lyrae stars (left panel of 
Fig.|9]l. Another way of looking at this is to assert that all the 
RR Lyrae variables in F2 come from M31, as stated above, so 
fiB = 14 ± \/l4. Then the number of observed RR Lyrae vari- 
ables in Fl is 17, so the number of RR Lyrae variables in Fl 
that do not come from M31 is ll-fj, B = 17- 14± v / 14 = 3 ±4 
(for integer numbers of stars). We therefore cannot claim to 
have detected RR Lyrae stars in M32 with reasonable confi- 
dence based on only Fl and F2, and we can put an upper limit 
of no more than 6 RR Lyrae belonging to M32 in Fl from this 
analysis. 

4.2. M32 RR Lyrae Population Inferred from F1-F4 

Unfortunately, both fields Fl and F2 suffer from small num- 
ber statistics due to the small angular size of the ACS/HRC. 
If we allow that Fl exceeds F2 just by chance, we should 
also ask if F2 exceeds the "true" background also by chance. 
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F435W-F555W F435W-F555W 

FIG. 4.— The (F435W-F555W, F555W) CMDs calibrated onto the HST VEGAMAG photometric system for fields Fl (left) and F2 (right). We show the 
location of the detected RR Lyrae variable stars. First-overtone (FO) and fundamental-mode (FU) pulsators are shown with empty and filled circles respectively. 
These CMDs are presented and discussed in M09. 
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FIG. 5. — As in Fig. [4] zoomed into the region of the detected RR Lyrae stars. Left: Fl. Right: F2. Numbers correspond to variables listed in Tables [2] (left 
panel) and|5](right panel), respectively. 



Fields F3 and F4 may provide some guidance here, as they 
cover much larger areas (as they were taken with ACS/WFC) 
and therefore have much smaller statistical errors, as shown in 
Figure [10] In addition these two fields nearly fall on the same 
M31 isophote that passes through Fl and F2. We therefore 
continue our search for M32 RR Lyrae variables assuming 
that fields F2-F4 constitute a fair sampling of the background. 

We must first consider the photometric and temporal com- 
pletess of F3 and F4. S09 did not perform ASTs on these 
fields to calculate the completeness, but accordin g to their lu- 
minos ity function and to their comparison with iBrown et al.l 
(2004), they assume 100% photometric completeness, which 
we will also assume here. Using the same procedure outlined 
in Section 3, we can estimate the temporal completeness for 
these fields. Figure QT| shows the computed probability for 
F3 and F4, revealing a slight difference between the sampling 
of the two fields F3 and F4. In particular for period values 
around (P a b) = 0.59 d, the probability decreases to 0.8 and 
0.95 when P = 0.65 and 0.55 d for F3 and F4, respectively. 



This is due to the shorter time window for these parallel fields 
than for the primary fields Fl and F2 (see Sec. 2.1). There 
are 45 (out of 324) RR Lyrae variables with periods 0.6-0.7 d 
(i.e., 80% temporal completeness) in F3 and 142 (out of 357) 
RR Lyrae variables with periods 0.5-0.6 d (i.e., 95% temporal 
completeness) in F4. We therefore estimate a total complete- 
ness of 97% and 98% for F3 and F4, respectively. 

Now we extend our equations ([T]!-© to include all three 
background fields, accounting for the different areas of the 
background fields and Fl . Because all of the background field 
measurements are independent, we can write the number of 
background counts as the sum of all three fields, corrected for 
completeness: 



F4 



i=F2 



A^ = £^ = 7L2 



(6) 



where is the number of stars in field i and c, is the com- 
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FIG. 6. — Atlas of light curves in the B and V bands of RR Lyrae stars detected in Fl. The B-band points have been shifted brighter by 2 mag for clarity. Error 
bars, as described in the text, take into account both the photometric errors as returned by DOLPHOT program as well as the scatter between the data and the 
model used to fit the Fourier series. The model is also shown in this Figure. For each variable star, its ID and period obtained by fitting the data points are shown 
in each panel as reported in columns 1 and 4 of Table|2] 



pleteness of field i. Then we can write equation (01 as 
P(H F \N B +N F ) oc J dfi B P(N B + N F \n F ,fi B ) 



PI '^b\N b > )P(fi F ), 



(7) 



where A B > /Ap\ is the ratio between the total area of the back- 
ground fields and the area of Fl, A B > /Apt = 2~2i=F2^B -Mfi = 
73. 

Assuming the same priors as in the previous section [equa- 
tion © and the following discussion], we find that the num- 
ber of RR Lyrae variables belonging to M32 in Fl has a most 
probable value of fi F = 7, with a 68% confidence interval of 
4-1 1 RR Lyrae variables (right panel of Fig. [9). Again, we 
can interpret this result as assuming that the estimated contri- 
bution from M31 RR Lyrae variables in Fl is 712/73 = 10, so 
the inferred number of M32 RR Lyrae variable in this field is 



17— 10 ± v 10 = 7 ± 3 (again, for integer numbers of stars). 
Note that this analysis is only correct if the surface bright- 
ness of M31 is constant across fields F1-F4. We can therefore 
claim to have detected seven RR Lyrae variables belonging to 
M32 in field Fl (with a ler upper limit of 1 1 RR Lyrae stars) 
under this assumption. Therefore we conclude that, indeed, 
M32 has an ancient population as represented by the detec- 
tion at 1-sigma level of RR Lyrae stars. 

4.3. Specific Frequency ofM32 RR Lyrae stars 

We can use this number of RR Lyrae stars belonging to 
M32 to estimate an upper limit on the specific frequency of 
M32 RR Lyrae stars (Srr) in Fl, assuming that its old, metal- 
poor 5 population resembles that of the Galactic globular clus- 
ters. Srr is defined as the number of RR Lyrae stars (Nrr) 
normalized to a total Galactic globular cluster luminosity of 

5 on the basis of the [Fe/H] value as derived in section |5. 3. II 
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FIG. 7. — As in Figure[6] but for RR Lyrae stars detected in F2. 



M Vt = -7.5 mag: 



'RR : 



„ r 1() (Mv, + 7.5)/2.5 



(8) 

(lHarrisI I1996I) . We therefore need to obtain the lumi- 
nosity of the old, metal-poor population of M32 in our 
field. The surface brightness of M32 in Fl is \iy w 
21.9 mag/D // (iKormendvet al.ll2009h . Assuming E(B - V) = 
0.08 {Schlegel et al. 1998) and a distance modulus (m-M)o = 
24.50 (M09) we obtain a total luminosity of M32 in Fl of 
My, ~ -10 mag. We then consider the metallicity distribu- 
tion function (MDF) of M32 derived by M09, from which a 
10 Gyr-old population having -2.3 < [Fe/H] < —1.3 consti- 
tutes 1 1 % of the total luminosity. The M32 luminosity of the 
old, metal-poor population in Fl is thus My = -7.6 mag. This 
implies Srr w 6.5 for this outer region of M32, with 68% 
confidence limits of 3.6 < Srr < 10, which is in reasonable 
agreement with Srr of Galactic globula r clusters with metal- 
licities [Fe/H] ~ -1.6 (see Fig. 10 in lBrown et alll2004l) . 

We can also estimate Srr for the stellar population of 
M31 sampled in F2. The surface brightness in this field is 



jtzy ~ 22.7mag/n", and assuming the same distance modu- 
lus and reddening than Fl, we obtain a total luminosity of 
M31 in F2 of My< ~ -9.24 mag. We consider again only the 
old, metal-poor population in F2 and calculate its luminosity. 
From the MDF of M3 1 by M09, a 10 Gyr-old population with 
metallicities -2.3 < [Fe/H] < —1.3 constitutes a 12% of the 
total luminosity, which translates into My, ~ -6.9 mag. The 
14 RR Lyrae stars found in F2 imply a Srr ~ 18 of M31 in 
this field. Note that this value is higher than Srr ~ 1 1 .2 esti- 
mated bv lBrown et all d2004t) for their M31 halo field. How- 
ever, the scatter in Srr in the Galactic globular cluster system 
is large enough to encompass these variations if M31's RR 
Lyrae population is simila r to the Milky Way's (cf. the discus- 
sion in iBrown et al.ll2004l) . 

One might be tempted to attempt to invert this analysis to 
determine the (upper limit on the) fraction of old, metal-poor 
stars in M32 and M3 1 . After all, we know the upper limit on 
the number of RR Lyrae stars we can associate with M32 in 
Fl and the number we can associate with M31 in F2. How- 
ever, no theoretical or empirical one-to-one relationship be- 
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TABLE 4 

Time-series magnitudes for RR Lyrae variables 



Julian date 
-2400000 



F435W 
mag 



B 

mag 



Julian date 
-2400000 



F555W 
mag 



Fl variable 1 



V 
mag 



53633.616 


25 


39±0 


07 


25 


44 


53635.480 


25 


42±0.07 


25 


48 


53633.631 


25 


22±0 


06 


25 


28 


53635.496 


25 


60±0.08 


25 


66 


53633.681 


25 


36±0 


07 


25 


42 


53635.545 


25 


72±0.08 


25 


78 


53633.697 


25 


69±0 


09 


25 


74 


53635.560 


25 


62±0.08 


25 


67 


53633.748 


25 


96±0 


10 


26 


01 


53635.612 


25 


44±0.06 


25 


50 


53633.764 


25 


89±0 


10 


25 


93 


53635.628 


25 


24±0.06 


25 


29 


53633.814 


25 


86±0 


10 


25 


90 


53635.678 


25 


25±0.06 


25 


31 


53633.830 


25 


50±0 


07 


25 


55 


53635.694 


25 


24±0.06 


25 


30 


53634.415 


25 


39±0 


08 


25 


44 


53636.213 


25 


17±0.06 


25 


22 


53634.430 


25 


37±0 


08 


25 


42 


53636.228 


25 


38±0.07 


25 


44 



Nrr 
80 
60 
40 



N 



RR 

C) 2 - 

80 - 
60 - 
40 - 
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- 


Fl 






- 






F2 










F3 
• 

1 , , 1 


F4 

i 








24 27 

A^M32 




30 



Fl 



F3 



F2 



F4 



NOTE. — The errors on HST VEGAMAG are the photometrical ones. The cal- 
ibration onto the Johnson-Cousins B- and V -bands as well as their errors have been 
discussed in the text. This is a sample table showing the format of the Table; the com- 
plete table can be found in the online version of the journal. 



£ 0.6 

5 
.a 

S 0.4 
a. 

0.2 



o.o L 
o.o 



1.0 

Period [days] 




Period [days] 



FIG. 8. — Spectra of the probability to detect variability in the period range 
from 0.2-2 d for the Fl (top) and F2 (bottom) fields observed with ACS/HRC 
(this paper). Note that the mean period found for Fl and F2 is 0.59 ± 0.1 1 d. 





FIG. 9. — Probability distribution functions for the expectation value of 
the number RR Lyrae variable stars belonging to M32 in field Fl. Left: 



P(N 



M32.FI 
RR Lyrae 



) = P(ij.f\Nb +Nf) [equation 0] using only the number of RR 

Lyrae stars found in Fl and F2. Right: PiN™ 2 ^) = P(p F \N B +N F ) [equa- 
tion (7)] using the number of RR Lyrae stars found in F1-F4, accounting for 
area differences and completeness. See text for more details. 
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FIG. 10. — The population of RR Lyrae stars in fields F1-F4 (see Fig. [T] 
for field placements). Top: RR Lyrae population (per square arcminute) as 
a function of the inferred M32 surface brightness (in the V band). Bottom: 
RR Lyrae population (per square arcminute) as a function of the projected 
distance from M32. 
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FIG. 1 1 . — Same as in Figure[8]but for the F3 and F4 fields observed with 
ACS/WFC and analyzed by S09. Periods found for F3 and F4 are 0.55 ± 
0.07 d and 0.56 ± 0.08 d, respectively. These spectra are not the same as in 
Figure [8] due to the fact that these fields were observed for only half of the 
time that Fl and F2 were. 

tween Srr and metallicity exists, due to the scatter in this 
SRR-[Fe/H] pl ane induced by the "second-para meter prob- 
lem" (see, e.g. JSuntzeff et al.lll99UlBrown et al.ll2004l) . Fur- 
thermore, we have used the known fraction of metal-poor stars 
in M32 in Fl (and in M31 in F2) determined by M09 and a 
reasonable guess for the ages of these metal-poor stars to esti- 
mate Srr, so such an analysis would be circular. Alternately, 
if there were a stellar evolution model that correctly predicted 
mass-loss along the first-ascent giant branch so that the distri- 
bution of stars along the horizontal branch was also correctly 
predicted, we could also perfor m this inversion; but no such 
model currently exists (see, e.g.. lSalaris & Cassisill2005h . 

4.4. Comparison with Alonso-Garcia et al. (2004) 
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TABLE 5 

Fields near M32 with claimed RR Lyrae detections 



ID 


RA 
(2000) 


DEC 
(2000) 


Time Window 




Instrument 


N(RR) M 


N(RR) TO 


N(RR)f 


FoV 
arcmin 2 


Completeness 


Fl a 


00:42:47.63 


40:50:27.4 


2005 Sep 20-24 (<- 


- 24 hr) 


ACS/HRC 


17 


13 


4 


0.25 


100% 


F2 a 


00:43:07.89 


40:54:14.5 


2006 Feb 6-12 (~. 


24 hr) 


ACS/HRC 


14 


8 


6 


0.25 


100% 


F3 b 


00:42:41.2 


40:46:38 


2005 Sep 22-24 (r. 


- 10 hr) 


ACS/WFC 


324 


267 


57 


9 


97% 


F4 b 


00:43:20.8 


40:57:25 


2006 Feb 9-12 (~ 


10 hr) 


ACSAVFC 


357 


288 


69 


9 


98% 


p 5 c 


00:43:01 


40:50:21 


1998 Nov 19 (~ 


4hr) 


WFPC2 


29 








5.7 


5-15% 


F6 C 


00:43:28 


41:03:14 


1998 Nov 20 (~ 


4hr) 


WFPC2 


16 








5.7 


5-15% 



d This paper 

b S09 

c lAIons o-Garci'a et al] 120041) 



> 25 




FIG. 12.— The JC-calibrated (V-I,V) CMD for field F5 and various com- 
pleteness levels (dashed lines). The arrows show the location where we ex- 
pect to find RR Lyrae stars. On this basis we conclude that the photometric 
completeness is between 20 and 75%. See M09 for more details. 



lAlonso-Ga rcfa et al] (120041) were the first to attempt to di- 
rectly detect an ancient population in M32 using RR Lyrae 
variables. They claimed to have detected 12 ± 8 RR Lyrae 
stars in a field further away from M32 than Fl. From this 
claimed detection, they conclude that 2.3% of M32's stellar 
population is old. 

In this section, we analyze the temporal and photometric 
completeness of Alonso-Garcfa et al.'s data in order to deter- 
mine if our detections and theirs are compatible. The authors 
do not give any estimate of the photometric completeness, but 
by examining their CMD we can see that the luminosity level 
of the horizontal branch is evidently at the detection limit of 
the WFPC2 exposures. We call their primary ("M32") field 
F5 and their background ("M31") field F6 in Table M09 
performed photometry and AST s on F5 to com pute their com- 
pleteness levels using HSTphot (Dolphin 2000). In Figure [T2l 
we show the {V-I, V) CMD calibrated onto t he JC photomet- 
ric sys tem using the calibrations provided by Holtzman et al. 
(1995). We conclude that the photometric completeness for 
RR Lyrae stars in F5 is 20-75%, depending strongly on color. 
Field F6 is similarly deep and therefore we assume that its 
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FIG. 13. — Same as in Figure[8]but for the F5 and F6 fields observed with 
WFPC2 and analyzed by Alonso-Gar ci'a et all <2004T) . 

photometric completeness is identical. Figure [13] shows that 
the temporal completeness computed for F5 and F6 using the 
method described in Section 3 is only as 0.2 for the mean pe- 
riod observed for RRab Lyrae (P ~ 0.59 days) and is 0.6- 
0.8 for RRc (P ~ 0.28-0.32 days). We stress here that the 
RRc detectability mostly depends on the quality of photom- 
etry due to the very low a mplitudes of FO pulsators and that 
lAlonso-Garcfa et all ( 120041) likely could not have detected any 
of these. Thus we can assume a total (photometric and tem- 
poral) completeness for these fields of around 5-15%. 

Given our own RR Lyrae detections in Fl, we can esti- 
mate their completeness in a different manner. We note that 
their WFPC2 fields are 22.8 times larger in area than our 
ACS/HRC fields, but that the M32 surface brightness in F5 
is /i^ 32 ' F5 ~ 23.7mag/D", about 5.2 times fainter than Fl. 
We assumed that the M31 surface brightness is the same in 
both fields, /i^ 131 ~ 22.7 mag/D". Then the total nu mber of 
RR Lyrae we predict that [Alonso-Garcfa et ail (12004ft should 
have detected if they had 100% completeness, given ATm3i.fi 
RR Lyrae variables belonging to M3 1 and ATm32,fi RR Lyrae 
variables belonging to M32 in Fl, is 



A/tot,F5 = 22.8ATm3i.fi + 22. 8ATm32.fi- 



jq-0.4/^(M32,F5) 



(9) 



I Q-0.4/j V '(M32,Fl) 

For ATm32.fi = 7^ RR Lyrae variables belonging to M32 in Fl, 
they should have detected AT to t,F5 = 258~gj RR Lyrae in total 
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(with fewer detected in F5 for more M32 variables in Fl). 
They claim to have detected 22 RR Lyrae variables in F5, so 
we estimate their complete ness to be cfs = 0.085^ jffit Due to 
this severe incompleteness. lAlonso-Garcfa et al] (120041) were 
not able to study the intrinsic properties of the stars they de- 
tected. In fact, they classified variable stars detected as RR 
Lyrae stars due only to their location in the CMD. 

We can now c ompute the number of M32 stars we should 
have seen in Fl if lAlonso-Garcfa et al.l (|2004) detected 12 ± 8 
RR Lyrae variables belonging to M32 in F5: 

A m l(T a4 ^ (FI) cpi 

AfcR,M32(Fl) - Nrr,M32(F5)— 10 _ 4w , (f5) — , (10) 

where Api and Ap5 are the areas of the two fields, /iy(Fl) 
and Hv(F5) are their surface brightnesses, and cpi and cp5 
are the completeness levels at the blue horizontal branch in 
the two fields. We assume here that this last ratio is 1 /0.085, 
based on the completeness estimate above. The scaling fac- 
tor in this equation is then 2.7, a number which is driven 
mostly by the incompleteness of F5. That is, we should 
have se en 2.7 times as many M3 2 RR Lyrae variables in 
Fl than lAlonso-Garcfa et al.l 00041 saw in F5, for a total of 
33 ± 22 RR Lyrae variables predicted to belong to M32 in 
our field. This is a factor of 4.6 higher than the number that 
we claim to have detected. However, the lower 68% confi- 
dence level of this estimate just overlaps with the upper 68% 
confidence level of our estimate of 7^ RR Lyrae variables 
predicted to belong to M32 in Fl. We can of course invert 
this procedure and predict how many M32 RR Lyrae variables 
lAlonso-Garcfa et all (120041) should have seen in F5 if our de- 
tection is correct. We predict that they should have seen 3^ 
RR Lyrae variables predicted to belong to M32 in F5, once 
again just in agreement (within the 68% confidence intervals) 
with their estim ate of 12 ±8. We therefore suggest that it is 
possible that Alons o-Garcfa et all ( 120041) detected bona fide 
M32 RR Lyrae variables in F5; but their severe incomplete- 
ness makes their detection highly uncertain. 

We summarize this section by stating that considering only 
our ACS/HRC fields Fl and F2, we have detected < 6 RR 
Lyrae variable stars belonging to M32 in Fl. Complementing 
our primary observations with our ACSAVFC parallel fields 
F3 and F4 allows a better determination of the likely M31 
background in Fl, and using this estimate, we claim to have 
detected 7^ RR Lyrae variable stars belonging to M32 in Fl, 
with a most probable value of seven variable stars. We there- 
fore conclude that M32 has an ancient population that can be 
detected without directly probing the oldest main-sequence 
turnoffs. 

5. RR LYRAE PROPERTIES 

We now ask whether we can separate the detected RR Lyrae 
stars belonging to the M32 and M3 1 stellar populations based 
on their intrinsic pulsation properties. 

5.1. RR Lyrae Star Colors and the Pulsation Instability Strip 

Does the average color of the RR Lyrae stars allow us to 
separate the two populations? In Figure [14] we show a blow- 
up of the CMDs in V calibrated onto the JC photometric sys- 
tem. Inspecting this figure we see a clear difference in the 
color ((B) - (V)) range of variables belonging to our different 
fields. In particular, Fl RR Lyrae colors are spread through a 
larger color range (0. 18 mag) than those in F2. 
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FIG. 14. — Zoomed view of the CMD around the positions of RR Lyrae 
stars. Mean magnitudes have been calibrated onto the Johnson-Cousin pho- 
tometric system. RR Lyrae stars from field Fl and those from F2 are dis- 
played as filled and empty circles, respectively. FO pulsators are emphasized 
by using other larger circles around the symbols used. The boundaries of the 
pulsation Instability Strip and the relative values we adopted to draw them 
are indicated in each panel, namely FOBEpi_p2 = 0.20 mag (left solid line), 
FREp2 = 0.40 mag (dashed line) and FREpi = 0.58 mag (right solid line). In- 
terestingly RR Lyrae stars in Fl appear redder than the ones in F2 suggesting 
a large spread in metallicity in this field (see the text for details). 

To understand this difference in color range, we need 
to discuss the mechanisms driving the radial stellar pulsa- 
tion. There are two main pulsation mechanisms, both re- 
lated to the opacity in the stellar envelope layers where 
partially-ionized elements are abundant (hydrogen and neu- 
tral or singly-ionized helium): the n- and 7- mechanisms 6 . 
These are directly correlated to the variations, in these lay- 
ers, of the opacity (k) and the adiabatic exponent (T3- l 7 ), 
respectively. These mechanisms can explain the pulsational 
instability strip (IS), the well-defined region in color where 
we find RR Lyrae stars and radial periodic pulsators in gen- 
eral (as classical and anomalous Cepheids, SX Phe variables, 
etc.) are observed. If a low-mass star (M < O.8M ) in the core 
helium- and hydrogen-shell-burning phase (the zero-age hor- 
izontal branch) crosses the IS during its evolution, it becomes 
a RR Lyrae variable star. 

Stars bluer than the First Overtone Blue Edge (FOBE) of 
the IS cannot pulsate because the ionization regions (located 
at almost constant temperatures) are too close to the surface, 
and thus the stellar envelope mass is too low to effectively 

6 These mechanisms are triggered in a region where an abundant ele- 
ment (hydrogen or helium) is partially ionized ( see Zhevakin 1953. [T^BH 
IBaker & Kippenhahr] ICoxlfF^l IBono & S tellingweri 1994"). Dur- 
ing an adiabatic c ompression the opacity increases with the temperature 
I King & Cox 1968), which is unlikely if the element is completely ionized. 
During this compression, heat is retained and the layer contributes to the 
instability of the structure. T his is called the K-mechanism and was first 
discussed by Eddington 1 1926), who called it the "valve mechanism." The 
7-mechanism allows, under an initial compression, the energy to be absorbed 
by the ionizing matter, instead of raising the local temperature, therefore low- 
ering the adiabatic exponent 7. The layer then tends to furthe r absorb heat 
during compression, leading to a driving force for the pulsations (Kin g & Coxl 
1968). It is evident that the two mechanisms are connected with each other. 
H, Hel and Hell layers are located at 13 000 K, 17,000 K, and 30,000-60,000 
K, respectively I Bono & Stellingweri 1994, and references therein), and thus 
this phenomenon involves only the stellar envelope. 

7 T 3 -l =<51og7y<51ogp 
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retain heat and act as a valve. Moving towards the red side, 
at lower effective temperatures, convection becomes more ef- 
ficient and, at the Funda mental Red Edge (FRE), quenches 
the pulsation mechanism (Smithl [l995h . The width of the IS is 
then just the color difference between the FOBE and the FRE. 
The FOBE location in the CMDs is well defined by theory, 
and we use this in Section 5.3.2 to find an ind ependent esti- 
mate of the distanc e modulus (as suggested by Caputo1 U997t 
Caput oet alj|200dl) . while the position of the FRE depends 
on the assumptions adopted to treat the convective transport, 
becoming bluer as the convective efficiency is increased. 

The IS depends on the initial stellar chemical composi- 
tion due to the dependence of the pulsation physics on the 
opacity and on the hydrogen and helium partial-ionization re- 
gions. In particular the metallicity and the convection effi- 
ciency (especially for the FRE) should have the largest effects 
on the boundaries of the IS. As the initial stellar metallic- 
ity increases, the opacity increases, resulting in a more effi- 
cient K-mechanism at lower effective temperature, which im- 
plies a redder FRE. The FOBE location is almost constant, 
as it depends essentially on the temperature of the region 
where helium is completely ionized (see, e.g.. lWalkerlll998t 
ICaputo et ai]l2000t [Fiorentino et afl l2002l) . 

In our case, shown in Figure [14] we do not see any differ- 
ence in the FOBE between the fields, with ((B) - (V))fobe = 
0.2 mag for both fields, whereas ((B) - (V))fre = 0.58 mag 
for Fl and 0.40 mag for F2. This evidence could be inter- 
preted as a suggestion that Fl RR Lyrae field stars may have 
a larger spread in metallicity or age than their M3 1 counter- 
parts, support for two different mixes of stellar populations in 
field Fl. Furthermore it is interestingly to note these red RR 
Lyrae stars (namely V6, V10, V14, V15, V16 and V17) in Fl 
seem to be closer to the M32 center, located as indicated by 
the arrow in Figure[2] Note also that the F2 RR Lyrae stars V3 
(the reddest RR e ) and V7 (the bluest RR fl i) define an empiri- 
cal "OR" region of the pulsation IS (see Figures|5]and in fl~4T> 
where, depending on their evolutionary path, RR Lyrae stars 
may pulsate in both FU or FO mode. 

5.1.1. RR Lyrae Reddening Evaluation 

To understand whether this difference in the color spreads is 
intrinsic — that is, due to the metallicities or ages of the stars — 
or just a re ddening effect, we m ust study the intrinsic redden- 
ing. From Schlegel et al. ( 1998) we know that Galactic fore- 
ground extinction in that direction is E(B-V) = 0.08 ±0.03 
mag and is essentially the same for both fields Fl and F2. 
However, we must also allow for extinction from the disk of 
M31, in case some of the RR Lyrae variables in Fl and F2 lie 
behind it. In the following we will use the intrinsic properties 
of the RR Lyrae stars to derive the total intrinsic reddening 
values for both fields in two independent ways. 

The first method is based on the insensitivity of ((B - 
V))fobe to the intrinsic properties of diff erent RR Lyra e 
populations, stressed for the first time by IWalkerl (119981) . 
Walker presents dereddened (magnitude-averaged) ((B - V)) 
IS boundary colors from accurate observations of nine Galac- 
tic and LMC globular clusters covering a range in metallicity 
from [Fe /H] = -1.1 to -2.2 dex. He concludes that the dered- 
dened color of the blue edge is at ((B-V))fobe = 0.18 ±0.02 
mag, with no discernible dependence on metallicity, while the 
color of the red edge shows a shift o f 0.04 ± 0.03 mag. This 
method was also used successfully by Cleme niini et al.l (120031) 
to estimate the average reddening value in LMC RR Lyrae 
stars. As discussed above both fields have the same FOBE 



color, and therefore we obtain the same estimate for their 
reddening. We find a reddened color of ((B) - (V))fobe 8 = 
0.20 ± 0.03 mag for both fields, as shown in Figure[l4] Thus, 
we infer a value of E(B-V) = 0.02 ± 0.04 mag for both fields. 
Interestingly, this value of the tota l reddening is more th an lcr 
smaller than the estimate from the lSchlegel et al.l(ll998l) map, 
which is meant to measure the foreground reddening. 

The second method was described originally by ISturchl 
( 1966). We stress here the importance of this method, as the 
calibration of the H I column dens ities with reddenin g is es- 
tablis hed on its basis and used by iBurstein & Heilesl (119781 
119821) . These authors ass umed an offse t of E(B-V) = 0.03 
from the values found by Sturch (1966). They fixed it un- 
der the assumpti on that at the Gal actic poles E(B-V) = 0, 
as suggested by iMcDonaldl (11977b in his analysis of RR 
Lyrae colors and H/3 indices, whereas Sturch found E(B- 
V) = 0.03. With this method we derive reddening for each 
RRab star from the (magnitude-averaged) color at minimum 
light (B-V) m i„ (phases between 0.5 and 0.7), the period P a i„ 
and the metal abundance [Fe/H] of the variables. The ap- 
plication of Sturch's method requires the knowledge of the 
metallicity of each individual RRab. Because we do not 
have this information from spectroscopy of these fields, we 
have decided to fix the average metallicity for both fi eld to 
the value [Fe/H] = -1.6 dex, as discussed in Section 15.3.11 
Sturch 's method was calibrated on Galactic field RR Lyrae 
stars (ISturchl 1966tlBlancol 19921) and has b een used on Galac- 
tic Globula r clusters dWalkeri|1990[ 119981) and on LMC RR 
Lyrae stars Clemen tini et al.l (120031) . returning, as expected, 
very good agreement with th e H 1-based redden i ng measure- 
ments taking into account the Burst ein & Heilesl (1 19781 [19821) 
extinction zero point. Walker ( 1992) used the following for- 
mulation of the Sturch's method, 

£(B-y) = (B-y), m -„-0.336-0.24P(d)-0.056[Fe/H], (H) 

where the reddening zero point has been adjusted to give 
E(B-V) = 0.0 mag at the Galactic poles, and the [Fe/H] is 
that of the Zinn & West (1984) metallicity scale. We in- 
fer mean reddening values of E(B-V) = 0.07 ±0.10 mag 
and 0.03 ± 0.08 mag for Fl and F2 resp ectively, in statistical 
agreement with each other and with the Schlegel et al.1 (119981) 
map. Note however that the scatter in the actual metallicities 
of the RR Lyraes will increase the scatter in these estimates. 

Depending on the assumed method, we therefore find 
two different estim ates of the reddening value: a) the 
iSchlegel et all (TT9981) map and the FOBE color method sug- 
gest the same reddening for both fields but different values de- 
pending on the method, 0.08 mag and 0.02 mag respectively; 
b) Sturch's method suggests a reddening slightly higher for 
Fl field than for F2 of about 0.04 mag. Taking into account 
the errors on both methods we do not find any significant dif- 
ference between the two evaluations as well as between the 
two fields. 

We will discuss all these cases when determining the dis- 
tance modulus: the firs t one assuming E(B-V ) = 0.08 for both 
fields by following the Schleg el et al.l (Q~998) map, the second 
one assuming E(B-V) = 0.02 mag for both fields as suggested 
by the FOBE method, and the third one assuming different 
values, e.g., E(B-V) = 0.07 mag for Fl and E(B-V) = 0.03 
mag for F2 found using Sturch's method. In all cases an ex- 

8 In our sample (Fl and F2) we verified that no significant difference is 
found between the FOBE as derived by the magnitude and intensity averaged 
colors, i.e. ((B)- <V»fobe-<(B-V)>fobe= 0.00 ±0.02 mag. 
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tinction law with Ry = 3.1 will be assumed. 

5.2. Mean Periods, Amplitudes and Oosterhoff Types 

In this section we focus our attention on their mean pe- 
riods, amplitudes and Oosterhoff types (Oo types) in order 
to distinguish between the two populations, M32 and M31 
halo (and/or outer disk) variables. Pulsational periods and 
amplitudes are of fundamental importance because they de- 
pend on the star structural parameters (mass, luminosity and 
effective temperature) and are distance- and reddening-free 
observables. 

RR Lyrae variables found in Galactic globular clusters can 
be divided in two distinct classes by the mean periods of 
the ir fundamental m ode pulsators: Oosterhoff type I and 
II dOosterhofj ri939f). In Oosterhoff type I (Ool) Galactic 
globular clusters, RRab variables have average periods of 
{P ah ) = 0.559 d, and in Oosterhoff type II (OoII) Galactic 
gl obular clusters they have average periods of (P a b) = 0.659 
d dClement et al. 2001), very few clusters have RRab with pe- 
riods in the range of 0.58 < (Pab) < 0.62 d (the "Oosterhoff 
gap"). These different Galactic Oosterhoff types may have 
resulted from different accretion and formation processes in 
the Galactic halo, a hypothesis supported by the difference 
in metallicities found for the two types (OoII RR Lyrae are on 
average more metal poor than Ool RR Lyrae, with ([Fe /H]) ss 
-1.8 and w -1.2, respectively, al though with large and over- 
lapping metallicity distributions: Szczy giel et alj|2009h . The 
Oosterhoff dichotomy may well hold a key to the formation 
history of the Galactic halo, but much evidence has suggested 
that this dichotomy cannot support the galaxy formation sce- 
nario in which present-day dwarf spheroidal (dSph) satellite 
galaxies of the Milky Way are the building blocks of our 
Galaxy (see, e.g.,ICatelan 2006). In fact the RR Lyrae vari- 
ables of these dSphs as well as those of their globular clus- 
ters fall prefere ntially into the "Oosterhoff gap." Furthermore, 
IClement et alJ (120011) investigated the properties of the RR 
Lyrae variables in Galactic globular clusters and found that 
the ratio between the number of RRc and RRab stars is con- 
stant for each of the two Oosterhoff types: N c /N lot!1 \ = 0.22 for 
the Ool clusters and N c /N totl ii = 0.48 for the OoII clusters. It 
is currently not clear whether the Oo dichotomy is a peculiar- 
ity of our own Galaxy; for example, the RR Lyrae stars ob- 
served in the M agellanic Clouds d o not appear to follow this 
di chotomy (seelAlcock et ail feOOO. and references therein). 

iBrown et al.l (120041) claimed that the RR Lyrae popula- 
tion they discovered in the M31 halo cannot be classified 
into either of the Oo types. They found a mean period of 
(P ah ) = 0.60 ±0.10 d (in the Oo gap) but N c /N lo ^ = 0.46, 
higher than the typical value for Ool clusters (0.22). On 
the other hand, Cle mentini et al.l (|2009) found that one of the 
most luminous M31 clusters, B514, shows a peculiarity in 
the Oosterhoff dichotomy. They found that the mean period 
of 82 RRab stars is (P ab ) = 0.58 d and N c /N tol!i i = 0.08, sug- 
gesting an Ool type, while at t he same time havi ng a very 
low metallicity, [Fe/H] ~ -1.8 dGalleti et al.ll2006l) . suggest- 
ing an OoII type. The v ery low value of the jV c /Motai ratio 
in this cluster found by Cle mentini et al.l d2009l) can be ex- 
plained by the incompleteness of their sample due to the in- 
trinsically low amplitudes of RRc stars. In a recent study 
based on parallel observations of our program taken with 
ACS/WFC, S09 found (Pab) = 0.55 ± 0.07 d and N c /N toa i = 
0.24 and (P ab ) = 0.56 ±0.08 d and N c /N tota i = 0.21 for the 
F3 and F4 fields (see Table [5), respectively, suggesting an 
Ool type classification. S09 applied the period-metallicity- 
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FIG. 15. — Bailey Diagram (V-band amplitude as a function of period) for 
RR Lyrae stars near M32 and in the disk of M31. RR Lyrae stars in Fl are 
shown in red, those in F2 are sho wn in blue, those found by S09 in F3 and F4 
are in grey, and those found in Brown et al. 1 2004) in F7 are in magenta. FU 
pulsators are circles and FO pulsators are triangles. As expected the FO pul- 
sators are concentrated in a region with low amplitude and periods, whereas 
the FU pulsators show a linear behavior, i.e., their amplitudes decrease as 
their periods increase. That trend is stressed by th e solid and das hed lines 
representing Oosterhoff types I and II, respectively I Clement 2000). Nearly 
all FU pulsators follow the Ool relation. 

amplitude relation to their large sample (almost 700 stars) and 
found TFe/Hl = -1.77 ± 0.06 for both fields, similar to the 
IClementini et al.l d2009l) results. 

As briefly discussed in Section [3] we find the same mean 
period (P a b) = 0.59 ±0.11 d and, within the uncertainties, 
the same ratios of FO to FU pulsators in the two fields: 
Nc/Ntou = 0.23^|] and N C /N M = 0.42^;| for Fl and F2, 
respectively. The error on the mean periods is just the stan- 
dard deviations of the average value computed on our very 
small sample of RRab Lyrae stars, i.e., 13 stars for Fl and 8 
for F2, and the errors on the ratio of FO to FU pulsators reflect 
only Poissonian counting uncertainties. 

We ask whether these properties are consistent with a sin- 
gle Oosterhoff ty pe. In Figur e [TBI w e show the composite 
Bailey diagram of IBrown et alJ d2004l) . S09, and our new de- 
tection of RR Lyrae variables. The comparison with other 
samples is made by assuming that their amplitudes, measured 
in the broad F606W filter (with a peak wavelength between 
Johnson-Cousins the V and R bands), underestima te the JC V- 
band amplitudes by i=s 8% (for further details see lBrown et alJ 
2004). Inspecting this figure, we see that RR Lyrae stars in 
both of our fields can be classified as Ool, as expected given 
the metallicity distribution of these fields (M09), and in par- 
ticular, their lack of well-developed blue horizontal branches. 

5.3. RR Lyrae Pulsation Relations 

We next attempt to disentangle the putative M32 RR Lyrae 
population from that of the M31 using various observa- 
tional manifestations of the theoretically-predicted RR Lyrae 
period-luminosity-temperature relation valid for every indi- 
vidual RR Lyrae star. 

5.3.1. RR Lyrae Metallicities: [Fe/H]-Period-Amplitude Relation 

To estimate the metallicity of t he RR Lyra e stars in the 
Large Magellanic Cloud (LMC), lAlcock et all ([2000) deter- 
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mined a relation between the metallicity [Fe/H] of an RRab 
variable, its period P„b, and its V-band amplitude Ay, 

[Fe/H] zty =-8.85[logP fli + 0.15A v ]-2.60, (12) 

where ZW refers to the lZinn & West! ([1984) metallicity scale. 
This relation has been calibrated from Galactic globular clus- 
ters with metallicities in the range of — 1 .4 < [Fe /H] zw < — 1. 1 
dex (and has been checked against the metallicities of Galactic 
field RR Lyraes). The uncertainty of this relation is Upe/H] = 
0.31. This relation shows very good agre ement with indepen - 
dent met allicity estimation s for the LMC (lAlcock et al.l200d|) . 
And IV dPritzl et alj l2002h . and And II dPritzl et aljhofrti) . 
Applying this relation to the RR Lyrae variables in our two 
fields we find mean values of ([Fe/H]pi) = — 1.52 ±0.10 ± 
0.48 dex and ([Fe/H] F2 ) = -1.65 ±0.11 ±0.40 dex. The first 
errors take into account the uncertainties of both periods and 
amplitudes, whereas the second errors are just the standard 
deviation from the mean value. There is clearly insufficient 
statistical evidence (given the estimated uncertainties) to con- 
clude that we are studying two completely different popula- 
tions, although the (insignificantly) higher metallicity of the 
RR Lyrae stars in Fl is in line with the results of M09. Note 
however that, as stated by lAlcock et al.l (120001) . this relation 
may be affected by evolutionary effects on the RRab variables 
and age effects in the stellar population. 

5.3.2. RR Lyrae Distance Modulus: Luminosity-Metallicity 
Relation and FOBE Method 

Much has been written about the luminosities of horizon- 
tal branch (HB) stars (particularly their V-band magnitudes) 
and their dependence on metallicity . Rece n t reviews of the 
subje ct inc lude papers bylChaboverl (119991) : ICacciaril (119991 
2003), and Cacciari & Clementini (2003). In this last paper, 
the authors describe the main properties of HB luminosities 
as follows: 
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• the relation M V (RR) = a [Fe/H] + (3 dSandagdll981allbl) 
can be considered linear, as a first approximation, with 
a variation of ~ 0.25 mag over 1 dex in metallicity; 

• at a given metallicity the evolution of low mass stars 
from the HB produces an intrinsic magnitude spread 
that can reach up to ~ 0.5 mag. This spread can in- 
crease as the metallic ity of the observed cluster in- 
creases dSandagd[T9 90): 

• the luminosity-metallicity relation is not strictly linear 
but depends on the HB morpholo gy and is related to 
the "second paramete r" problem (Ca puto et al.1 l2000l 
iDemarque etaTll2000l) . 

RR Lyrae variables remain however excellent distance in- 
dicators once these effects are properly known and taken into 
account. We assume a linear My (RR)-[Fe/H] relation, the av- 
erage of several methods: 

M V {RR) = (0.23 ±0.04)[Fe/H] + (0.93 ±0.12) (13) 

(for further details see ICacciari & Clementinil [20031) . To de- 
rive the distance modulus, we used this equation with ob- 
served mean magnitudes of (V) =25.34±0.15magforFl and 
(V) = 25.30±0.12 mag forF2 and a metallicity of [Fe/H] = 
-1.6 dex for both fields (roughly the mean of our inferred 
metallicities in Sec. 5.3.1). Furthermore, we assumed differ- 
ent values for th e reddening values as derived and discussed 
in Section IB". 1.1 1 
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FIG. 16. — RR Lyrae variable stars in the (V)-logP plane, where the visu- 
alization of the FOBE method is shown for an independent dist ance modulus 
(DM) determination. The symbols are the same as in Figure [T4l The solid line 
represents the FOBE edge, by assuming a stellar mass value of M = 0.7Mq 
and a metallicity [Fe/H] = —1.6. With these assumptions we find a DM of 
fio = 24.56 for a reddening value of 0.08 mag for both fields. If we assume 
the minimu m redd ening estimate of 0.02 mag as we derived (FOBE method, 
see Section I5.1.1K all the stars move towards fainter magnitudes an amount 
indicated by the arrow. In this case, we obtain a DM of fio = 24.72 mag 
(shown by the dashed line). The dotted lines represent the FOBE edges we 
obtained by taking into account a 0. 1 Mq spread in stellar mass around the 
initial assumed value. 

(a) ^o(Fl) = 24.53 ±0.21 mag and /i (F2) = 24.49 ±0.19 mag 

by us ing E(B-V) = 0.08 m ag for both fields, according 
to the ISchlegel et all dl998l) map; 

(b) jU (Fl) = 24.56 ±0.21 mag and /i (F2) = 24.64 ±0.19 mag 

by using E(B-V) = 0.07 mag and E(B -V) = 03 re- 
spectively for Fl and F2, according to the lSturchl (jl966) 
method; and 

(c) ^o(Fl) = 24.72±0.21magandAt (F2) = 24.59±0.19mag 

by using E(B-V) = 0.02 mag for both fields, according 
to the FOBE method. 

The errors are conservative and take into account the scatter 
in the mean (V) value, metallicity uncertainties (of about 0. 10 
dex), and the intrinsic uncertainty in the Mv(RR)-[Fe/H] re- 
lation. Note that only taking into account for a differential 
reddening between the two fields (case b) can one find Fl (the 
field closer to M32) to be slightly closer to us than F2. As 
shown by S09, mostly because of the large distance of M3 1 
from us and the relatively short distance between M31 and 
M32, these distance moduli do not provide any further infor- 
mation to disentangle the two populations. These moduli ap- 
pear to be in very go od agreement with estimates of both M3 1 
(un = 24.44 ± 0.1 1. iFreedman & Madorelll990l: u = 24.5 ± 



O.l.lBrown et al.l2004t un = 24.47±0.07 
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; 24.46 ±0.11, 



120051: mo = 24.54 ± 0.07. lSaha et aIll2006T^ _ 
S09) and M3 2 (un = 24.55 ± 0.0 8, ifonrv et al.l l200lt /z = 
24.39 ±0.08, Uensen et al] 120031) distance moduli and with 
a recent estimate of the distance to M32 obtained with our 
data using the Red Clump method, /io(Fl) = 24.50 ± 0. 12 mag 
(M09). 

The FOBE met hod described in ICaputol d 19971) and 
Caput oet al.l d2000) provides an independent distance estima- 
tor. It is a graphical method based on the predicted period- 
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FIG. 17. — RR Lyrae stars in the reddening-free Wesenhei t pla ne: (V) — 
3.1 X (B — V) vs. logP. The symbols are the same as in Figure [l4l Solid and 
dashed lines represent the linear fits to Fl and F2 datasets respectively. The 
slopes of these relationships show a small but insignificant difference, as the 
Fl relationship is well within the lcr and 3cr uncertainties of the F2 relation 
(thin dashed and dotted lines, respectively). 

luminosity (PL) relation for pulsators located along the first- 
overtone blue edge (FOBE) and seems quite robust for clus- 
ters with significant numbers of RRc variables. Using this 
procedure, one obtains a distribution of the cluster RR Lyrae 
stars in the My-logP plane once a distance modulus has been 
assumed. By matching the observed distribution of RRc vari- 
ables with the follo wing theoretical relation for the FOBE 
dCaputo et al.ll2000l) . 

My (FOBE) = -0.685(±0.027) - 2.255 logP(FOBE) 

-1.2591og(M/M ) + O.O581ogZ, (14) 

we can obtain an independent estimate of the distance mod- 
ulus. We assume the mean metallicities derived above, 
[Fe/H] = —1.6 for both fields, corresponding to logZ = -3.3, 
and we take M = 0.7M© from evolutionary HB models 
for RRc variables , with an uncertainty of the order of 4% 
dBono et al.ll2.Q03b . The FOBE method then yields a distance 
modulus of /io = 24.56 ±0.10 mag for both fields assuming 
E(B-V) = 0.08 mag. If we also take into acc ount th e min- 
imum reddening value we obtained in Section IB. 1.1 1 of 0.02 
mag, the distance modulus increases to /io = 24.72 ±0.10 
mag. 

5.3.3. PLC and PLC -Amplitude Relations 

Our final attempt to disentangle the two populations is to 
consider at the same time all the information on the RR Lyrae 
variables we have obtained (including magnitudes, ampli- 
tudes, colors and periods) that reflect the individual properties 
of these stars. In Figures[l7]and[T8]we show our detected RR 
Lyrae stars in the reddening-free Wesenheit plane, (V) -Ry(< 
B > - <V >)-jAy vs period, without (7 = 0) or with (7 ^ 0) 
taking into account a dependence on the stellar amplitudes, re- 
spectively. Because our sample of RR Lyrae variables is very 
small, we "fundamentalized" the first-overtone (FO) pulsators 
using the relation log Pf = logfVo + 0.13, where log Pf is the 
fundamentalized period ( Ivan Albada & Bakerlll973l) . In Fig- 
ure [17] we see that both RR Lyrae groups are located in the 
same region; the solid and dashed lines represent the best lin- 
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FIG. 18. — RR Lyrae stars in the reddening-free Wesenheit plane, taking 
into account the V-band amplitudes: (V) — 3 .1 X (< B > - < V >) — jAy vs 
log P. The symbols are the same as in Figure [T"4l Solid and dashed solid lines 
represent the fits to our data obtained with +O.OI7 and —0.297 values for Fl 
and F2, respectively. There is again an insignificant difference between the 
two fits, as shown in the lcr and 3cr uncertainties for the F2 relation (thin 
dashed and dotted lines, respectively). 

ear fit across the data points for Fl and F2, respectively. As 
we can see only the slope of the relation changes, but this 
slope difference is statistically insignificant. The thin dashed 
and dotted lines represent the dispersion around the relation 
for F2, for la and 3a respectively. In Figure [18] we see that 
the periods of RR Lyrae stars in F2 do not show any depen- 
dence on their own amplitudes, whereas the periods of RR 
Lyrae stars in Fl show a mild dependence on their pulsation 
amplitudes of 7 = -0.19 dex. However, the two relationships 
are in statistical agreement with each other, suggesting that 
the RR Lyrae properties of the two groups are very similar. 

We summarize this section by stating that the properties of 
RR Lyrae variables in our two fields, Fl and F2, are statis- 
tically inseparable. The only difference we find between the 
two fields is in the color of the FRE, where a population of 
Fl RRab variables are redder than the F2 RRab variables. We 
suggest that these stars may in fact belong to the M32 RR 
Lyrae population, but the small number of detected RR Lyrae 
stars prevent us from making a definitive claim as to the true 
cause of this color difference. 

6. WHAT HAVE WE LEARNED ABOUT M32, AND WHAT DO WE 
STILL NOT KNOW? 

In this paper we present ACS/HRC observations of fields 
near M32 (and supplemental ACS/WFC observations) to 
search for RR Lyrae variable stars. The detection of RR Lyrae 
variable stars represents the only way to obtain information 
about the presence of an ancient, metal -poor population (> 10 
Gyr and [Fe/H] < — 1 dex) in M32 from optical data not deep 
enough to detect the oldest MSTO stars. 

We have detected 17 RR Lyrae variable stars in Fl, our 
field closest to M32, and 14 RR Lyrae variable stars in F2, 
our background field in M3 1 's disk. We can only claim to 
have detected an upper limit of 6 RR Lyrae stars belonging to 
M32 in Fl based on these two fields alone. We can better con- 
strain the M32 RR Lyrae population by extending our analysis 
to our ACS/WFC parallel fields F3 and F4 and by assuming 
that the M31 surface brightness is constant over these fields. 
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With this extension, we can claim to have detected 7^ RR 
Lyrae variable stars belonging to M32 in Fl, and therefore 
we claim that M32 has an ancient population. The implied 
specific frequency of RR Lyrae stars in M32 is Srr « 6.5, 
with a 68% confidence interval of 3.6 < Srr < 10, similar to 
typical intermediate-metallicity ([Fe/H] ~ —1.5 dex) Galac- 
tic globular clusters. In making these estimates, note that we 
have assumed that the metal-poor population constitutes 11% 
of the total luminosity of M32 in Fl following the metallic - 
ity distribution function presented in M09 and further that this 
population has an age of 10 Gyr. 

We have also used the background M3 1 population to in- 
fer a specific frequency of RR Lyrae variables in M31's disk, 
Srr » 18, that is higher than the value found bv lBrown et alj 
(2004) in M31's halo, Srr k, 11. This may reflect merely 
small fluctuations in the horizontal-branch morphology with 
position in M3 1 or it may reflect some property of M3 1 's disk; 
it is difficult to state conclusively the cause without a better 
understanding of how the horizontal branch is populated (and 
therefore how mass is lost on the first-ascent giant branch). 

Even though we claim to have detected bona fide M32 
RR Lyrae variable stars in Fl, the pulsational properties 
stars in fields Fl and F2 present nearly no significant dif- 
ferences. They have indistinguishable mean V magnitudes 
[(V) = 25.34 ±0.15 mag and /i = 24.53 ± 0.21 mag forFl; 
(V) = 25.30±0.12 mag and /i = 24.49±0.19 mag for F2 
by assuming E(B-V)=0.08 mag and [Fe/H] = —1.6 dex for 
both fields], the same mean periods ((P a b) = 0.59 ±0.11 d), 
and the same distribution in the Bailey (V-band amplitude-log 
Period) diagram, and insignificantly different ratios of RRc to 
RRab types (N c /N ab = 0.30^ and N c /N ab = 0.75^ respec- 
tively). By using the relation between the metallicity [Fe/H] 
of each RRab variable, its period P ab , and its V-band ampli- 
tude Ay, we find mean values of ([Fe/H])pi =-1.52±0.10 
dex and ([Fe/H]) F2 = -1.65 ±0.11 dex. We stress here that 



these metallicities, as well as those found by S09 ([Fe/H] = 
-1.77 ± 0.06 from nearly 700 stars), are far away from the 
"solar-like" metallicity used by iBrown et alJ (l2000h to inter- 
pret the UV excess they found in the core of M32 as a hot 
horizontal branch. However, we do find a population of Fl 
RRab variables that are redder than F2 RRab variables. We 
suggest that these stars may in fact belong to the M32 RR 
Lyrae population, but the small-number statistics imposed by 
the small ACS/HRC area makes a definitive claim difficult or 
even impossible. 

We have detected M32's long sought-after ancient, metal- 
poor population through its RR Lyrae population. But we find 
that this RR Lyrae population is nearly indistinguishable in 
its mean pulsational properties from M31, and for that mat- 
ter, its specific frequency. Does this imply that M32 and M3 1 
formed, evolved, and (likely) interacted in such a way that 
their ancient, metal-poor populations share some commonal- 
ity? Or is this just a coincidence? We hope that further ex- 
ploration of M32's RR Lyrae population, from fields closer 
to the center of M32, and high-resolution spectroscopy of in- 
dividual stars, revealing their detailed chemical compositions 
and therefore details about their formation histories, in M32 
and the disk of M3 1 will eventually become available to help 
answer these questions. 
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